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Abstract 

Battery  life  is  an  important,  yet  technically  challenging,  issue  for  battery  development  and  application.  Adequately  estimating  battery  life 
requires  a  significant  amount  of  testing  and  modeling  effort  to  validate  the  results.  Integrated  battery  testing  and  modeling  is  quite  feasible 
today  to  simulate  battery  performance,  and  therefore  applicable  to  predict  its  life.  A  relatively  simple  equivalent-circuit  model  (ECM)  is  used 
in  this  work  to  show  that  such  an  integrated  approach  can  actually  lead  to  a  high-fidelity  simulation  of  a  lithium-ion  cell’s  performance  and 
life.  The  methodology  to  model  the  cell’s  capacity  fade  during  thermal  aging  is  described  to  illustrate  its  applicability  to  battery  calendar  life 
prediction. 
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1.  Introduction 

Accurate  prediction  of  a  battery’s  life,  either  the  calendar 
or  the  cycle  life,  is  a  great  technical  challenge  to  any  battery 
application.  Battery  life  is  an  important  issue  for  both  trac¬ 
tion  and  stationary  applications  because  it  is  critically  related 
to  battery  reliability  and  dependability,  which  in  turn  deter¬ 
mines  a  power  source’s  quality  and  eventual  life-cycle  cost. 
A  battery’s  actual  service  life,  however,  depends  on  its  his¬ 
tory,  during  both  the  storage/standby  and  mission/duty  pe¬ 
riods  experienced  by  the  battery  through  its  lifetime.  Any 
sensible  approach  to  predict  battery  life  therefore  has  to  ad¬ 
dress  impacts  from  both  the  storage/standby  and  mission/duty 
periods. 

In  the  past  decade,  due  to  substantial  improvements  in 
computer  computation  power  and  software  capability,  bat¬ 
tery  modeling  and  simulation  [1-8]  has  enjoyed  significant 
advancements.  At  the  same  time,  experimental  techniques 
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that  allow  detailed  investigation  of  interfacial  and  bulk  prop¬ 
erties  of  electrode  materials  have  contributed  to  a  better  un¬ 
derstanding  of  cell  performance  and  degradation  [9] .  It  is  fea¬ 
sible  to  develop  an  integrated  battery  testing  and  simulation 
capability  to  assist  battery  R&D  and  operation  [5,6]. 

There  are  a  few  attempts  in  the  past  to  predict  lithium- 
ion  battery  capacity.  For  example,  Fuller  et  al.  [10]  used  a 
‘first-principles’  electrochemical  model  to  estimate  lithium- 
polymer  cell  capacity.  Rakhmatov  et  al.  [11]  proposed  an 
analytical  model  for  lithium-ion  cells  used  in  portable  elec¬ 
tronic  systems  that  can  predict  battery  lifetime.  Spotnitz  [12] 
incorporated  SEI  growth  into  Fuller’s  model  and  began  to 
look  into  the  correlation  of  impedance  change  with  capacity 
fade.  Ramadass  et  al.  [13]  attempted  to  incorporate  solvent 
reduction  reaction  into  their  first-principles  electrochemical 
model  to  predict  capacity  fade. 

We  are  much  interested  in  taking  a  comprehensive,  concur¬ 
rent  and  hybrid  approach  [7,8]  to  develop  tools  and  strategies 
that  can  predict  the  end-of-life  (EOL)  of  a  battery.  Determina¬ 
tion  of  the  EOL  of  a  battery  system  in  service  is  very  difficult, 
expensive,  and  often  destructive  to  the  system.  Therefore, 
non-invasive,  non-destructive  techniques  that  can  determine 
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Fig.  1.  The  ECM  used  in  this  work. 


the  EOL  without  removing  the  system  from  service  are  highly 
desired. 

Sandia  has  been  involved  in  studying  the  performance  of 
a  group  of  18650-size  lithium-ion  batteries  (LIB)  that  was 
fabricated  for  the  US  Department  of  Energy  (USDOE)  Ad¬ 
vanced  Technology  Development  ( ATD)  Program.  These  bat¬ 
teries  contain  a  high-power  chemistry  [9,14-16].  One  of  the 
objectives  of  the  ATD  Program  is  to  develop  procedures  to 
make  rapid  comparisons  of  performance-degradation  rates 
and  predictions  of  battery  life.  This  paper  will  discuss  a  sim¬ 
ple  LIB  modeling  approach  that  uses  an  equivalent-circuit 
model  (ECM)  to  simulate  cell  performance,  particularly  the 
capacity  fade  phenomenon.  This  example  shows  that  we  can 
simulate  battery  performance  changes  due  to  thermal  aging, 
which  is  one  of  the  most  influential  factors  impacting  battery 
calendar  life  during  storage,  standby  or  operation  periods. 


2.  The  equivalent-circuit  model  (ECM) 

There  are  several  ECM  approaches  of  different  nature  and 
flavor  reported  by  others  (e.g.,  [17-19])  for  simulating  LIB 
performance.  A  schematic  of  the  ECM  used  in  this  work  is 
shown  in  Fig.  1.  The  model  resembles  to  that  used  by  Ver- 
brugge  and  Conell  [19]  for  Ni-MH  cells.  We  favor  this  model 
due  to  its  simplicity,  yet  flexibility,  in  describing  an  electro¬ 
chemical  system  via  the  separation  of  all  ohmic  resistance 
components  from  all  faradic  non-linear  components,  as  they 
are  lumped  into  R\  and  R^C,  respectively.  This  generic  na¬ 
ture  makes  this  ECM  applicable  to  a  variety  of  chemistries 
(e.g.,  valve-regulated  lead  acid  [7],  Ni-MH  [19],  or  LIB  in 
this  paper). 

To  construct  a  valid  model,  we  first  need  to  incorporate  the 
state-of-charge  (SOC)-dependent  open  circuit  voltage  (OCV) 
and  resistance  ( R\  and  R2)  values  into  the  model.  Fig.  2 
shows  these  SOC-dependent  values  for  the  chemistry  used 
in  this  work.  The  OCV  values  were  determined  from  a  cell 
discharged  at  C/25  rate,  as  reported  in  [20].  The  resistance 
versus  SOC  relationship  can  be  derived  from  two  sources  of 
experimental  data,  either  from  ac  impedance  measurements 
or  from  the  galvanostatic  dc  polarization,  which  yields  dis¬ 
charge  curves  at  different  rates.  A  more  detailed  discussion  of 
how  the  values  were  derived  can  be  found  in  [8].  We  should 
note  that  in  the  model,  the  resistance  R2  was  assumed  to  con¬ 
sist  of  two  independent  contributions,  R 2  and  7?",  as  shown  in 
Fig.  2,  to  approximate  the  resistance  change  with  SOC.  The 
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Fig.  2.  SOC-dependent  OCV  and  resistance  values  used  in  ECM. 

function  of  R2  follows  a  power  law  with  SOC,  while  that  of 
R2  is  an  exponential  function  with  SOC,  as  shown  in  Eqs.  (1) 
and  (2). 

R'2  =  a  +  b(  SOCf  (1) 

R2  =  dexp[(l  -SOC)<?]  (2) 

Thus, 

R2  =  a  +  b( SOCf  +  d  exp[(l  -  SOC)e]  (3) 

These  functions  are  empirical  and  do  not  imply  any  physi¬ 
cal  meaning  in  this  work,  although  further  analysis  of  their 
dependence  with  aging  time  and  conditions  could  provide  us 
some  clues  as  to  what  process  might  be  shaping  the  physical 
change  of  the  cell  performance.  This  would  need  to  be  ver¬ 
ified  by  additional  experimental  evidence.  The  capacitance 
value,  C,  can  be  derived  from  the  characteristic  frequency, 
fC2  =  l/R2C,  which  is  yielded  from  the  Nyquist  plot  of  the 
cell  as  demonstrated  in  Fig.  3.  The  time-dependent  cell  volt¬ 
age  discharge  curve  can  then  be  calculated  from  the  ECM  for 
a  constant-current  condition  [19],  according  to  Eq.  (4): 

V(t)  =  9/p.Q-t/RiC  +  _  IRl  _  /fl2(i  _  e-'/fi2c)  (4) 


Fig.  3.  Nyquist  plot  for  cell  #319  measured  at  60%  and  100%  SOC. 
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where  I  is  constant,  g(0)  is  the  nominal  capacity,  and  V0  the 
nominal  SOC-dependent  cell  OCV.  Thus,  for  any  given  time 
step,  the  charge  passed  is  calculated  to  derive  the  SOC  value 
at  the  end  of  the  step,  and  the  associated  OCV  and  resistance 
values  are  plugged  into  Eq.  (4)  to  give  the  cell  voltage  as  a 
function  of  time. 

3.  Experimental 

Quallion  LLC  (Sylmar,  CA)  provided  test  cells  as  part  of 
the  USDOE  ATD  Program  efforts,  and  Sandia  along  with  four 
other  National  Labs  was  commissioned  to  evaluate  the  cell 
performance  for  hybrid  electric  vehicle  applications.  Detailed 
descriptions  of  the  chemistry,  cell  configuration,  test  proto¬ 
cols  and  procedures,  and  test  results  can  be  found  in  [9, 1 6,20] . 
In  short,  this  high-power  cell  chemistry  uses  a  cathode  con¬ 
sisting  of  LiNi0.8Co0.15Al0.05O2,  an  anode  fabricated  with 
MAG- 10  graphite,  and  an  electrolyte  consisting  of  ethylene 
carbonate  (EC)/ethyl  methyl  carbonate  (EMC)  (3:7  wt.%  ra¬ 
tio)  with  1.2  M  LiPF6. 

The  data  used  in  this  paper  were  extracted  from  that  ob¬ 
tained  for  a  cell  (#319)  in  the  pool  of  Sandia’s  test  cells 
as  an  example.  This  particular  cell  went  through  thermal 
aging  at  55  °C  and  100%  SOC  for  about  20  weeks,  while 
the  other  lithium-ion  cells  were  aged  at  temperatures  rang¬ 
ing  from  25  °C  to  55  °C  and  at  60%,  80%,  or  100%  SOC 
to  examine  the  cell’s  power  and  capacity  degradation.  Cell 
#319  was  chosen  because  its  behavior  is  representative  of 
the  group  of  cells  that  went  through  the  same  thermal-aging 
conditions. 

The  cell  performance  as  received  and  during  subsequent 
tests  at  four- week  intervals  was  evaluated  by  a  reference  per¬ 
formance  test  (RPT),  as  described  in  [20].  The  cell  capacity 
was  determined  at  C/1  and  C/25  discharge  rates  in  RPTs  after 
each  thermal-aging  period,  ac  impedance  data  were  also  col¬ 
lected,  and  were  intended  to  identify  which  cell  components 
were  changing  in  the  aging  process.  Other  National  Labora¬ 
tories  on  the  ATD  Program  are  conducting  more  detailed  di¬ 
agnostic  analyses  to  understand  the  degradation  mechanism. 
The  discharge  curves  from  the  discharge  regimes  at  C/25  and 
Cl  1  and  the  ac  complex  impedance  data  generated  in  the  RPT s 
were  used  in  deriving  the  parameters  for  model  development 
and  validation. 


4.  Results  and  discussion 

Fig.  4  shows  a  series  of  discharge  curves  at  various  rates, 
from  C/25  to  10C,  simulated  by  the  ECM  using  Eq.  (4).  Also 
shown  are  two  sets  of  experimental  data  obtained  at  C/25 
and  C/1  rate.  The  very  high  degree  of  agreement  between  the 
actual  data  and  the  simulated  results  gives  us  a  high  level  of 
confidence  about  the  resistance  values  that  we  used  in  the 
model.  It  is  worth  mentioning  that  the  fit  of  the  parameters 
is  also  demonstrated  by  the  ability  to  capture  the  essence  of 
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Fig.  4.  Discharge  curves  simulated  by  ECM  using  Eq.  (4)  and  the  data  ob¬ 
tained  at  C/25  and  C/1  for  comparison  and  validation. 

the  shape  of  the  discharge  curves.  This  is  achieved  primarily 
by  the  use  of  the  two  contributions,  R'?  and  R ",  in  R2 .  The 
more  consistency  achieved,  the  better  we  can  distinguish  any 
subtle  differences  from  the  parameters  used  to  interpret  the 
results.  Eventually,  this  capability  can  be  used  to  analyze  sub¬ 
tle  differences  from  various  phenomenological  assumptions 
and  the  corresponding  reflection  in  the  parameters  applied  to 
the  battery’s  physical  model  to  help  us  study  the  impacts  on 
the  power-delivery  capability,  capacity,  and  life. 

The  ability  to  simulate  cell  voltage  under  galvanostatic 
charge  or  discharge  regimes  is  quite  useful  for  battery  R&D. 
For  example,  by  imposing  a  suitable  cut-off  condition,  we 
can  calculate  the  amount  of  charge  put  in  and  released  from 
the  cell  to  yield  the  rate  capacity,  as  well  as  charge  efficiency. 
Furthermore,  if  the  aging  effects  and  degradation  rates  are 
known,  we  can  perform  simulations  using  various  tempera¬ 
ture  and  power  conditions  imposed  onto  the  cell  to  simulate 
the  life  performance. 

In  order  to  simulate  the  capacity  fade  of  the  cell,  we  need 
to  know  how  cell  characteristics  such  as  OCV  and  resistance 
values  change  with  aging  time  and  conditions.  Since  OCV 
and  R\  are  intrinsic  to  the  cell,  their  changes  under  thermal 
aging  should  be  negligible.  Therefore,  we  treated  them  as 
invariant  in  this  work.  R\  =  0.01 1  £2  was  used  throughout  the 
work,  independent  of  SOC  and  aging  conditions. 

The  only  characteristic  that  requires  attention  is  the  change 
in  R2.  Fig.  5  shows  how  the  five  parameters  in  Eq.  (3)  change 
with  the  aging  period  under  the  same  aging  condition.  These 
relationships  allow  us  to  interpolate  or  extrapolate  parame¬ 
ter  values  over  a  specific  aging  period  to  permit  prediction 
of  the  cell  performance  characteristics  under  the  aging  con¬ 
dition.  Parameter  a  is  a  constant  and  independent  of  SOC, 
and  parameters  b  and  c  change  in  linear  fashions.  Parameter 
d ,  the  pre-exponential  factor,  changes  in  a  3rd  order  polyno¬ 
mial  while  e  changes  in  a  power  law  fashion. 

Fig.  6  shows  total  cell  resistance  from  the  above  five  pa¬ 
rameters  in  Eq.  (3)  and  R\.  The  overall  resistance  of  the  cell 
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Fig.  5.  Trends  of  parameters  a  to  e  in  Eq.  (3)  versus  aging  period. 


as  a  function  of  SOC  increases  as  the  cell  ages.  We  also  no¬ 
tice  that  although  the  resistance  increases  occurred  across  the 
entire  SOC  range,  they  are  more  noticeable  in  the  lower  SOC 
region,  thus  leading  to  a  higher  rate  of  capacity  fade.  The 
separation  of  R2  into  two  independent  contributions  of  R2 
and  R 2  allows  a  disproportional  rate  of  increase  in  R'?  and 
R2,  making  the  approach  more  suitable  and  easier  to  reflect 
characteristic  changes  in  capacity  fade. 

Fig.  7  shows  the  results  of  using  this  set  of  parameters 
in  ECM  to  simulate  the  discharge  behavior  of  the  cell  going 
through  the  thermal-aging  experiments.  The  five  discharge 
curves  simulated  from  the  ECM  are  consistent  with  the  exper¬ 
imental  data.  The  differences  between  the  predicted  capacity 
and  the  experimental  data  are  compared  in  Fig.  8. 
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Fig.  6.  Characteristics  of  resistance  change  versus  SOC  as  the  cell  goes 
through  the  thermal  aging.  Symbols  are  experimental  data.  Lines  are  values 
estimated  from  Eq.  (3)  with  parameters  in  Fig.  5  and  R\ . 


SOC,  % 


Fig.  7.  Simulated  discharge  curves  (symbols)  with  two  independent  con¬ 
tributions  in  R2  in  the  ECM,  in  comparison  with  the  experimental  results 
(lines). 
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Fig.  8.  Comparison  of  capacity  fade  measured  from  experiments  and  sim¬ 
ulated  from  ECM.  Polynomial  fitting  allows  projection  of  battery  calendar 
life  under  the  thermal-aging  condition. 


B.Y.  Liaw  et  al.  /  Journal  of  Power  Sources  140  (2005)  157-161 


161 


The  introduction  of  the  two  resistance  contributions  in  R2 , 
one  that  follows  a  power  law  dependence  with  SOC,  and  the 
other  exponential  is  somehow  essential  in  order  to  properly 
simulate  the  discharge  behavior,  and  thus  the  capacity  change 
induced  by  the  thermal-aging  process.  Although  not  specif¬ 
ically  illustrated,  we  have  also  tried  to  treat  the  resistance 
R2  with  a  single  expression  of  polynomial  fitting.  This  treat¬ 
ment,  however,  would  not  lead  to  a  disproportional  increase 
in  resistance  at  different  regions  of  the  SOC.  Thus,  the  pro¬ 
found  capacity  fade  in  the  lower  SOC  region  could  not  be 
captured  with  this  particular  approach  using  a  consistent  set 
of  parameterization.  In  other  words,  unlike  the  clear  relation¬ 
ships  displayed  in  Fig.  5,  we  might,  at  best,  generate  each 
curve  under  degradation  with  an  isolated  set  of  purely  em¬ 
pirical  parameters,  such  that  it  bears  no  analytical  meaning 
at  all.  Therefore,  the  introduction  of  the  two  resistance  con¬ 
tributions  in  R2  not  only  allows  us  to  assign  an  asymmetric 
dependence  of  R2  with  SOC  but  also  enables  us  to  retain  a 
consistent  trend  of  parameterization  for  capacity  fade  predic¬ 
tion. 

From  the  parameterization  shown  in  Fig.  5,  we  would  fur¬ 
ther  like  to  point  out  that 

(i)  The  a  +  b( SOC) c  term  in  Eq.  ( 1 )  did  change  with  thermal 
aging  and  induced  an  increase  in  polarization  resistance 
that  undermined  the  cell  capacity.  However,  this  increase 
only  attributed  to  the  capacity  loss  in  the  high  SOC  re¬ 
gion.  It  did  not  affect  the  capacity  fade  in  a  predominant 
manner  since  the  capacity  fade  primarily  came  from  the 
capacity  loss  in  the  low  SOC  region. 

(ii)  The  capacity  fade  is  primarily  governed  by  the  term 
<iexp[(l— SOC)p]  in  Eq.  (2),  and  profoundly  following 
the  rapid  increase  in  the  pre-exponential  factor,  d. 

(iii)  Finally,  Fig.  8  shows  the  model  prediction  and  actual 
experimental  data,  in  which  the  capacity  fade  projec¬ 
tion  and  measurement  are  compared.  The  polynomial 
fittings  will  allow  the  prediction  of  battery  calendar  life 
under  this  particular  thermal-aging  process,  once  the 
EOL  condition  is  defined. 


5.  Conclusion 

We  have  demonstrated  that  high-fidelity  simulation  of  bat¬ 
tery  performance  and  capacity  fade  can  be  achieved  with  a 
simple  equivalent-circuit  model  using  a  consistent  set  of  pa¬ 
rameters  that  reflect  thermal  aging.  This  high-performance 
simulation  can  assist  us  to  predict  the  battery  life  under 
thermal-aging  conditions,  which  are  the  most  critical  fac¬ 
tors  in  determining  the  battery  calendar  life.  With  increas¬ 
ing  accuracy  in  predicting  the  battery  performance  and 
degradation,  we  can  use  this  approach  to  understand  the 
trend  of  resistance  change  during  thermal  aging  and  the 
impact  from  such  a  change  on  battery  performance  and 
life. 
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